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1.  Research  Objectives 


A  deep  iindersr;indlng  of  viscoelasticity  and  rheology  is  crucial  to  advanced  materials 
engineering  and  process  design.  Examples  of  such  advanced  materials  are  high-strength 
polymers  and  additi\es  for  lubricants;  process  design  problems  include  spinning  of  s3-n- 
thetic  fibers  and  injection  molding.  The  materials  involved  in  these  technologies  are  often 
highly  elastic  and  very  viscous.  As  a  consecpience.  they  often  disphij'  behavior  in'^erme- 
diate  between  that  of  a  solid  and  that  of  a  fluid,  and  their  djmamic  response  involves 
multiple  time-scales.  The  understanding  of  the  equations  of  motion  coupled  with  various 
constitutive  assumptions  at  the  mathematical  level  is  crucial  for  modelling,  design  of  al¬ 
gorithms  and  computational  solution  of  practical  problems.  The  goals  of  our  ri^search  can 
be  stiminarized  as  follows; 

1.  To  understand  qualitative  properties  of  the  mathematical  model;  the  global  existence 
and  uniqueness  of  solutions,  dependence  on  data,  regularitj"  and  asvmiptotic  behavior 
of  solutions  for  large  time,  approach  to  s^'eady  states,  etc. 

2.  To  understand  the  pht'sical  consec[uences  of  the  model.  Do  am'  or  all  of  the  possible 
solutions  make  phr'sical  sense?  Do  solutions  that  have  mathematically  interesting 
ch.ara,ctcr  coiTespoud  to  observed  phenomena?  Do  they  predict  l^ehavior  that  should 
be  studied  in  the  laboratory'?  What  solutions  to  the  problem  are  relevant  to  processing 
and  design? 

3.  To  understand  the  physical  model.  How  do  the  observed  solutions  correspond  to  the 
molecular  or  continuum  model  on  which  they  are  based?  Can  the  character  of  these 
solutions  serve  to  validate  the  physical  model  or  suggest  improvements  in  it? 

4.  To  design  numerical  methods  that  account  for  the  mathematics  and  reproduce  the 
physics.  The  current  focus  is  on  one-dimensir*".-'i  -nodel  problems,  but  in  the  long¬ 
term.  the  goal  is  to  adapt  results  for  the  mode,  "'•r  ’  Icms  to  flows  in  comple.x  geome¬ 
tries. 

5.  To  study  the  broad  mathematical  implications  of  1.  -  4.  for  these  and  related  classes 
of  problems.  We  note  that  concepts  such  as  weak  solutions,  shocks,  shear  bands,  phase 
changes,  etc  arise  in  a  variety  of  problems  of  physical  and  technological  importance. 


2.  Research  Accomplishments  of  Principal  Investigators 

In  recent  joint  work  with  B.  .1.  Plohr  [10  -  13.  24,  25].  we  investigated  several  models 
wuh  tkc.-c  features.  We  discovered  several  dramatic  new  phenomena  in  shear  flow  that 
appear  to  be  relevant  to  material  instabilities  that  disrupt  polymer  processes. 

One  striking  phenomenon  has  been  oljserved  by  Vinogradov  et  ai.  (1972*)  in  the  flow 
of  monodisperse  polyisoprenes  through  capillaries.  They  found  that  the  volumetric  flow¬ 
rate  increased  dramatically  at  a  critical  stress  that  w^as  independent  of  molecular  weight. 

-r.T-iookvdi  oi  dismissed  by  luvologists 


This  phenomenon,  ’.'.duck  ij  callcc 


because  no  plausible  mechanism  was  known  to  explain  it  in  the  context  of  steady  flow's. 
Spurt  was  lumped  together  with  instabilities  such  as  ‘‘slip.”  “apparent  slip.”  and  “melt 
fracture.”  which  are  poorly  understood  (see  recent  surveys  by  M.  Denn  (19S9)  and  .1. 
Pearson  (1985)).  While  regarded  as  anomalous,  these  instabilities  can  severeh'  disrupt 


*  References  identified  Iw  date  are  listed  in  chronological  order  in  Sec.  5 
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polynii?!'  procesbes:  they  can  be  avoided  in  practice  only  with  ad  hoc  engineering  expedients. 
The  mechanisms  of  such  phenomena  were  not  understood  because  the  governing  equations 
are  analytically  intractable,  and  because  popular  numerical  methods  for  steady  viscoelastic 
fluid  flows  falter  in  this  “high  Weissenberg  number”  regime,  arid  thus  cannot  model  the 
spurt  phenomenon. 

By  contrast,  our  recent  work  shows  that  satisfactory  modelling  and  explanation  of 
spurt  and  related  phenomena  require  studying  the  full  dynamics  of  the  equations  of  motion 
and  constitutii-e  equations.  The  common  and  key  feature  of  constitutive  models  that 
exhibit  spurt  and  related  phenomena  is  a  non-monotonic  relation  between  the  steady  shear 
stress  and  strain  rate.  This  allows  jumps  in  the  steady  strain  rate  to  form  when  the  driving 
pressure  gradient  exceeds  a  critical  value;  such  jumps  correspond  to  the  sudden  increase 
in  volumetric  flow  rate  observed  in  the  experiments  of  Vinogradov,  et  al.  (1972).  The 
governing  systems  used  to  model  such  one- dimensional  flows  are  analyzed  in  [6.  10  -  12]  by 
numerical  techniques  and  simulation,  and  in  our  more  recent  work  by  analytical  methods. 
The  systems  derive  from  fully  three-dimensional  differential  constitutive  relations  with  m- 
relrxation  times  (based  on  work  of  Johnson  and  Segalman  (1977)  and  Oldroyd  (1958)). 
They  are  evolutionary,  globally  well  posed,  and  they  possess  discontinuous  steady  states 
of  the  type  mentioned  above  that  lead  to  an  explanation  of  new  phenomena. 

We  model  these  flows  by  decomposing  the  total  shear  stress  into  a  polymer  contribu¬ 
tion.  evolving  in  accordance  with  a  differential  constitutive  relation  with  a  single  relaxation 
time  and  a  Newtonian  viscosity  contribution.  The  flows  can  also  be  modelled  by  a  system 
based  on  a  differential  constitutive  law  with  two  widely  spaced  relaxation  times  but  no 
Newtonian  viscosity  contribution.  We  developed  three  numerical  methods  for  tlie  model 
initial  boundary  value  problem;  one  based  on  Plohr’s  recent  work  on  Riemann  problems 
for  hyperbolic  conservation  laws  with  a  regularizing  Newtonian  viscosity,  another  based 
on  using  the  total  stress  (rather  than  the  velocity  gradient)  as  a  dependent  variable  in  the 
balance  of  linear  momentum,  and  a  tliird  based  on  appropriate  modifications  of  an  exist¬ 
ing  solid-mechanics  code.  Each  of  the  three  techniques  yields  satisfactory  qualitative  and 
quantitative  agreement  with  Vinogradov’s  experimental  results.  Numerical  simulation  [6. 
10  -  12]  of  transient  flows  at  high  Weissenberg  (Deborah)  number  and  very  low  Reynolds 
number  exhibited  spurt,  shape  memory,  and  hysteresis;  furthermore,  it  predicted  other 
effects,  such  as  latency,  normal  stress  oscillations,  and  molecular  weight  dependence  of 
hysteresis,  that  should  be  analysed  further  and  tested  in  rheological  experiment. 

Our  numerical  work  suggested  a  scaling  valid  for  Vinogradov’s  material  which  allowed 
the  governing  system  of  partial  differential  equations  to  be  reduced  to  a  pair  of  nonlinear 
ordin-^ry  differential  equations  by  setting  a  key  parameter  (ratio  of  Reynolds  number  to 
Deborah  number  which  is  of  order  10“^^  for  a  Vinogradov  material)  to  zero.  The  resulting 
(Ivmunics  can  be  ilerermined  completely,  i  h^*  results  of  the  phase  plane  analysis  are 
summarized  in  [13];  analytical  results  will  appear  in  [24,  25).  .4.t  the  mathematical  level,  it 
remains  a  challenge  (being  addressed)  to  resolve  the  singular  perturbation  problem  arising 
in  this  situation  between  the  full  and  reduced  systems.  In  [25],  we  also  justify  the  use  of 
Newtonian  viscosity  to  model  the  effects  of  relaxation  times  which  are  very  much  shorter 
than  the  fundiuaental  time  constant  of  the  macroinolecules.  We  determine  conditions  under 
which  the  character  of  the  phase  portrait  of  a  more  realistic,  two  relaxation-time  model 
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without  Newtonian  viscosity  matches  matches  the  phase  portrait  of  the  simpler  model  with 
Newtonian  viscosity.  Our  analytic  results  helped  us  to  understand  certain  rather  dramatic 
effects  that  were  observed  in  the  numerical  solutions.  Moreover,  they  enabled  us  to  identify 
the  parametric  dependence  of  these  effects  and  to  predict  their  experimental  signature. 

The  analytical  results  discussed  above  ha,ve  motivated  Nohel  (in  collaboration  with 
A.  Tzavaras  and  R.  Pego  )  to  undertake  a  deeper  analytical  study  of  the  full  unsteady 
shear-flow  problem  keeping  the  inertial  term.  In  [26.  27].  we  discuss  results  on  the  nonlinear 
stability  of  discontinuous  steady  states  of  a  model  initial-boundary  value  problem  in  one 
space  dimension  for  incompressible,  isothermal  shear  flow  of  a  non-Newtonian  fluid  between 
parallel  plates,  driven  by  a  constant  pressure  gradient.  The  model  incorporates  several 
properties  of  the  more  complex  problem  described  above.  The  non-Newtonian  contribution 
to  the  shear  stress  is  assumed  to  satisfy  a  simple  differential  constitutive  law.  Again,  the 
key  feature  is  a  non-monotone  relation  between  the  total  steady  shear  stress  and  steady 
shear  strain  rate  that  results  in  steady  states  having,  in  general,  discontinuities  in  the  strain 
rate.  We  e.xplain  why  every  solution  tends  to  a  steady  state  as  t  — >  oo.  and  we  identify 
steady  states  that  are  stable.  In  current  work.  Nohel,  Pego  and  Tzavaras  attempting  to 
extend  these  results  to  the  full  problem  discussed  above.  Preliminary  research  indicates 
that  discontinuous  steady  states  of  the  more  complex  problem  are  stable  (perhaps  not 
asymptotically)  if  the  ratio  of  Reynolds  number  to  Deborah  number  is  sufficiently  small. 

To  achieve  our  goals,  our  interdisciplinary  program  has  adapted  and  extended  tools  in 
nonlinear  analysis  of  partial  differential  equations,  analytical  and  computational  techniques 
for  hyperbolic  conservation  laws,  and  computational  techniques  from  nonlinear  structural 
dyniimics.  .Awareness  of  the  latest  developments  in  the  rheology  of  viscoelastic  liquids  was 
vital  in  these  efforts. 

3.  Research  Accomplishments  of  Associates  Students 

(A)  A.  Tzavaras,  Van  Vleck  Assistant  Professor.  Shear  instabilities  in  the  form  of  shear 
bands  are  often  observed  when  metals  are  deformed  at  high  strain  rates.  Since  shear  bands 
diminish  the  strength  of  materials  and  are  often  precursors  to  rupture,  their  understanding 
is  critical  for  the  development  of  improved  materials.  According  to  a  popular  theory,  shear 
band  formation  is  the  result  of  a  destabilizing  feedback  mechanism  induced  by  thermal 
softening  properties  of  materials.  .At  high  strain  rates,  non-tmiform  straining  induces  non- 
uniform  heating  which,  in  turn,  enhances  the  plastic  flow  at  hotter  regions  and  reduces 
it  at  colder  regions.  This  creates  a  destabilizing  feedback  mechanism  to  winch  there  is 
opposition  from  internal  dissipation  and  from  strain  hardening  properties  of  the  material. 
In  order  to  assess  the  contribution  of  the  above  factors  we  have  been  studying  the  adiabatic 
plastic  shearing  of  an  infinite  plate  subjected  to  prescribed  tractions  or  steady  shearing  at 
the  boundmies.  For  constitutive  law  we  choose  a  power  law  with  parameters  measuring 
the  relative  weight  of  thermal  softening,  strain  hardenining  and  strain  rate  sensitivity. 
The  analysis  of  the  above  problem  leads  to  the  study  of  a  system  of  partial  differential 
equations  consisting  of  a  parabolic  equation  coupled  through  the  diffusion  coefficient  with 
two  equations  of  hyperbolic  type.  The  main  question  is  to  study  whether  the  solution 
stabilizes  at  large  times  or  whether  nonuniformities  develop  and  the  material  exhibits 
an  unstable  response.  This  question  is  studied  in  [34]  with  the  following  results.  In 


case  the  sheaiing  deformation  is  caused  by  prescribed  tractions  the  parameter  space  can 
be  decomposed  into  two  subregions.  In  one  of  them  the  solutions  are  asymptotically 
attracted  to  a  "uniform”  dynamic  solution,  while  in  the  complement  instabilities  develop. 
The  analysis  indicates  that  the  latter  are  associated  with  a  collapse  of  the  ability  of  the 
material  to  diriuse  the  applied  stress. 

In  a  joint  program  with  M.  Slemrod  we  are  studying  the  wave  fan  admissibility  crite¬ 
rion  for  the  solution  of  the  Rieinaiin  problem  for  hyperbolic  conservation  laws.  It  is  well 
known  that  in  solving  the  Riemann  problem  one  encounters  a  loss  of  uniqueness  that  has 
to  be  accounted  for  by  choosing  m  apimopriate  admissibility  criterion.  According  to  the 
wave  fan  admissibility  criterion  the  admissible  solutions  are  chosen  as  limits  of  appropriate 
"viscosity  regailarized”  problems,  that  are  rigged  so  as  to  preserve  the  invariance  of  the 
equations  under  rescaling  of  the  independent  variables.  This  approach  has  been  succesfully 
tested  for  the  Riemann  problem  for  the  equations  of  gas  dynamics  in  one  space  dimension 
in  Eulerian  coordinates  [32], 

(B)  Deborah  Brandon,  Post-Doctoral  Research  Associate.  The  research  involves  the 
modeling  and  analysis  of  problems  in  heat  conduction  in  one  dimensional  materials  with 
memory.  During  the  reporting  period,  work  was  completed  on  the  construction  of  concrete 
integral  models  for  such  materials,  based  on  a  very  general  theory  developed  by  Gurtin  and 
Pipkin  which  is  compatible  with  the  second  law  of  thermodynamics  but  which  leaves  the 
form  of  the  constitutive  functionals  for  the  heat  flux  and  for  the  internal  energ\-  unspecified. 
The  results  will  appear  in  [3].  This  theory  which  ignores  mechanistic  effects  predicts  a 
finite  speed  of  propagation  for  thermal  disturbances.  However,  Gurtin  and  Pipkin  did 
ncjt  provide  any  concrete  examples.  Her  principal  achievement  is  the  explicit  construction 
of  constitutive  functionals  for  the  heat  flux  and  internal  energy  which  depend  on  athe 
temporal  history  of  the  temperature  gradient  and  possibly  on  the  present  value  as  well 
as  on  the  history  of  the  temperature  and  which  also  satisfy  the  neccssarj'  and  sufficient 
conditions  of  the  Gurtin-Pipkin  theory.  These  functionals  are  motivated  by  models  used  in 
viscoelasticity.  The  key  difficulty  is  that  the  crucial  condition  which  must  be  understood 
deeply  and  then  properly  interpreted  is  in  the  form  of  a  functional  differential  (Clausius- 
Duhem)  inequality.  As  a  consequence,  she  shows  that  compatibility  with  the  second  law 
requires  the  interned  energy  (as  well  as  the  heat  flux)  to  have  a  nontrivial  dependence  on 
the  histor}'  of  the  temperature  gradient  as  well  as  on  the  present  value  and  the  history  of 
temperature. 

(C)  R.  C.  Rogers,  Van  Vleck  Assistant  Professor.  In  addition  to  the  joint  work  with 
Nohel  and  Tzavaras  [15]  described  above,  Rogers  persued  research  on  self-  effect  problems 
in  electro-magneto-elasticity  [17,  21],  modelling  ferromagnetic  materials  [IS],  and  problems 
in  compensated  compactness  and  weak  continuity  [19,  20], 

(D)  K.  Barki,  Ph.D.  student.  Mr.  Barki’s  work  [35]  is  part  of  a  project  involving 
Malkus.  Barki.  Y.-C.  Tsai,  and  another  student.  R.  Cornwell  (who  is  supported  by  AFOSR 
funds  from  a  different  grant.  .4FOSR  89-0220).  The  pvirpose  of  the  project  is  to  develop 
a  dynamic  finite  element  code  for  the  .solution  of  problems  in  solid  thermo- vdscoplasticity 
and  the  use  of  this  code  to  investigate  the  calibration  and  resulting  resolution  of  laboratory 
stress-analysis  devices  based  on  differential  thermography.  Malkus  is  collaborating  with 
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M.  Plesha.  R.  Rowlands,  and  B.  Saiidor  (Dept.  Engineering  Mechs..  U.-W.)  on  a  new 
device,  called  SPATE,  that  promises  to  be  a  powerful  tool  in  expeiimental  stress  analysis. 
The  device  operates  by  sensing  and  interpreting  small  temperature  differences  induced  by 
c}'clic  loading.  At  present,  the  device  has  several  operating  restrictions;  first,  the  theory 
behind  it  applies  only  to  linearly  elastic  materials;  second,  the  output  is  the  trace  of  the 
stress  tensor,  not  individual  stress  components;  and  third,  this  stress  invariant  can  be 
sensed  reliably  only  at  the  surface  of  the  test  specimen.  To  make  SPATE  a  more  flexible 
tool  in  stress  analysis,  it  is  important  (i)  to  determine  principal  stress  components  from 
the  data  for  the  stress  invariant  and  (ii)  to  use  the  separated  surface  stress  distributions  to 
deduce  interior  stress  distributions.  Solving  these  problems  requires  an  accurate  model  of 
material  behavior  and  appropriate  boundary  conditions,  so  that  the  governing  equations 
may  be  formulated  and  solved. 

Malkus  and  Cornwell  have  successfully  simulated  thermograplric  signals  emitted  by  a 
viscoelastic  body  in  uniaxial,  cyclic  loading  by  developing  algorithms  for  coupled  thermo¬ 
viscoelasticity,  including  dissipation  heating.  Extensions  to  include  materials  with  yield 
stresses  are  in  progress.  We  are  now  developing  the  capability  to  solve  thermally-dependent 
problems  in  a  higher  dimensional  setting,  based  on  the  work  of  K.  Barki.  His  thesis  re¬ 
search  involves  the  development  of  new  transient  algorithms  for  finite  element  analysis  of 
viscoplastic  materials.  The  major  focus  is  on  a  class  of  implicit /explicit  methods  with 
origins  in  structural  dynamics,  generalizing  the  method  developed  by  Malkus  to  solve  dy¬ 
namic  spurt.  There  are  two  fundamental  ideas:  first  to  allow  implicit  treatment  of  the 
elastic  (highest)  wave  speed,  which  avoids  severe  stability  limits;  second,  to  avoid  the  need 
for  stress  interpolations  by  advancing  the  constitutive  equations  at  Gauss-points  in  ODE 
fashion.  Mr.  Barki  is  focu.ssing  on  problems  in  solid.  3-D  viscophtsticity  in  a  Lagrangian 
frame,  but  the  ideas  seem  to  be  useful  for  fluid  flow  problems  in  an  Eulerian  frame.  This 
idea  is  being  pursued  by  Malkus.  Mr.  Barki  has  been  using  the  CMS  VAX  system  and 
the  SDSC  Cray  XMP/48.  He  expects  his  Ph.D.  in  the  fall  of  19S9. 

(E)  Y.-C,  Tsai,  Ph.D.  Candidate.  Mr.  Tsai  has  been  investigating  the  theoretical  me¬ 
chanics  of  thermally  coupled  viscoplasticity  and  the  experimental  determination  of  model 
parameters;  of  particular  interest  is  the  relaxation  spectrum  of  the  heat  capacity  term  in 
the  energy  equation.  The  next  phase  of  the  planned  research  is  to  investigate  the  con¬ 
sequence  of  viscoplastic  behavior  on  differential  thermography,  by  simulating  the  SPATE 
stress  analysis  system.  The  theory  behind  dilferential  thermography  only  rigorously  applies 
to  linearly  elastic  material,  but  material  failure  mechanisms  often  involve  plastic  deforma¬ 
tion.  stress  relaxation,  and  creep.  How  such  mechanisms  would  be  sensed  by  differential 
thermograph}':  we  are  addressing  the  question  of  how  the  device  could  be  calibrated,  and 
what  the  accuracy  of  such  measurements  would  be. 

4.  Publications 

A.  BOOKS 

1.  R.  Cook.  D.  S.  Malkus.  and  M.  E.  Plesha.  Concepts  and  .4pplications  of  Finite  Element 

.4nalysi.s.  Third  Edition,  John  Wiley  and  Sons.  New  York.  19S9.  630  pp. 

2.  M.  Renardy.  W.J.  Hrusa.  and  J..4.  Nohel,  Mathematical  Problems  in  Viscoelasticity, 

Pitman  Monographs  and  Surveys  in  Pure  and  Applied  Mathematics,  vol.  35,  Longman 
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Scientific  Sc  Technical.  Essex,  England  (10S7),  273  pp 

B.  PAPERS  PUBLISHED  OR  ACCEPTED 

3.  D.  Brandon  and  W.  J.  Hrusa,  Construction  of  a  class  of  integral  models  for  heat  flow 
in  materials  with  memory,  Jonmai  of  Integral  Equations,  accepted. 

4.  A.  Friedman  and  A.  E.  Tzavaras,  Combustion  in  a  porous  medium,  SIAM  J.  Math. 
Anal.  19  (19SS),  509-519. 

5.  W.  J.  Hrusa.  J,  A.  Nohel  and  M.  Renardy.  Initial  value  problems  in  viscoelasticity. 
Applied  Mechanics  Reviews  41  (19SS),  371-37S. 

6.  R.  W.  Kolkka,  D.  S.  Malkus,  M.  G.  Hansen.  G.  R.  lerley.  and  R.  A.  Worthing.  Spurt 
phenomena  of  the  Johnson  Alodel  fluid  and  related  models,  J.  Non-Newtonian  Fluid 
Mechanics  29  (19SS).  303-325. 

7.  D.  S.  Malkus  and  M.  F.  Webster.  On  the  accuracy  of  finite  element  and  finite  difference 
predictions  of  non-Newtonian  slot  pressures  for  a  Maxwell  fluid.  J.  Non-Newtonian 
Fluid  Mechanics  25  (19S7),  93-127. 

S.  D.  S.  Malkus  and  X.  Qiu,  Divisor  structure  of  finite  element  eigenproblems  arising 
from  negative  and  zero  masses.  Coinput.  Meths.  Appl.  Mech.  Engrg.  66  (19SS). 
365-36S. 

9.  D  S.  Malkus,  M.  E.  Plesha,  and  M.-R.  Liu.  Reversed  stabifity  conditions  in  transient 
finite  element  analysis.  Comput.  Meths.  Appl.  Mech.  Engrg.  68  (19SS).  97-114. 

10.  D.  S.  Malkus,  J.  A.  Nohel,  and  B.  J.  Plohr,  Time-dependent  shear  flow  of  a  non- 
Newtonian  fluid  (preliminary  report),  Transactions  of  the  Sixth  Army  Conference  on 
Applied  Mathematics  and  Computing  (Boulder.  1088). 

11.  D.  S.  Malk'is.  J.  A.  Nohel,  and  B.  J.  Plohr,  Tiine-dependent  shear  flow  of  a  non- 
Newtcnian  fluid.  Current  Problems  in  Hyberbolic  Problems:  Riemann  Problems  and 
Computations  (Bowdoin,1988),  Contemporary  Mathematics.  Amer.  Math.  Soc.,  ed. 
B.  Lindquist,  accepted. 

12.  D.  S.  Malkus,  J.  A.  Nohel.  and  B.  J.  Plohr.  Dynamics  of  shear  flow  of  a  non-Newtonian 
fluid.  Journal  of  Computational  Physics,  accepted. 

13.  D.  S.  Malkus.  J.  A,  Nohel.  and  B.  J.  Plohr.  Anah  ■sis  of  spurt  phenomena  for  a  non- 
Newtonian  fluid.  Proceedings  of  Conference  on  Problems  that  Change  Type  (Stutt¬ 
gart.  October  (1088).  Springer  Verlag  Lecture  Notes,  accepted. 

14.  J.A.  Nchel  and  M.  Renardy,  Development  of  singularities  in  nonlinear  viscoelasticity, 
Proceedings  of  Workshop  on  Amorphous  Polymers,  IMA  Volumes  in  Mathematics  and 
its  Applications.  Vol.  6,  Springer- Verlag  (1987),  139-152. 

15.  J.A.  Nohel.  R.C.  Rogers,  and  A.  Tzavai'as.  Weak  solutions  for  a  nonlinear  system  in 
viscoelasticity.  Comm,  in  P.D.E.  13  (1988),  97-127. 

16.  J.A.  Nohel.  R.C.  Rogers,  and  A.E.  Tzavaras.  Hyperbolic  conservation  laws  in  vis¬ 
coelasticity.  in:  Volterra  Integro  Differential  Equations  in  Banach  Spaces  and  Appli¬ 
cations,  G.  da  Prato  S:  M.  lanelli  (Eds),  Longman  Scientific  and  Technical  (in  press). 

17.  Robert  C.  Rogers  and  Stuart  S.  Antman.  Steady-state  problems  of  nonlinear  electro¬ 
magneto-thermo-elasticity.  Arch.  Rat.  Mech.  Anal.  95  (1986),  279-323. 

IS.  Robert  C.  Rogers,  Nonlocal  problems  in  electromagnetism,  in  Stuart  S.  Antman,  J.L. 
Ericksen.  David  Kinderlehrer,  and  Ingo  Muller,  editors.  Metastability  and  Incom¬ 
pletely  Posed  Problems,  IMA,  Springer- Verlag,  New  York  (1986). 
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19.  Joel  \V.  Robbia.  Robert  C.  Rogers,  and  Blake  Temple.  On  weak  continuity  and  the 
Hodge  decomposition.  Tnui-sa ctions  of  the  AMS  302  (10S7). 

20.  Robert  C.  Rogers  and  Blake  Tennple.  A  sufficient  condition  for  weak  continuity  of 
polynoniicil.s  in  the  method  of  compensated  compactness,  Transnctions  of  the  AMS. 
accepted. 

21.  Robert  C.  Rogers.  Nonlocal  variational  i^roblems  in  nonlinear  electro-magneto-elas- 
ticit  v.  SIAM  J.  of  Math.  Analysis,  accepted. 

22.  .4.  E.  Tzavaras.  Effect  of  thermal  softening  in  shearing  of  strain  rate  dependent  ma¬ 
terials.  Arch.  Rational  Mech.  and  Anal.  90  (19S7).  349-374. 

23.  M.  ’^'ao  and  D.  S.  Malkus.  Boundary  node  correction  and  super-convergence  in  the 
finite  element  method,  preprint  (19S9).  Int.  J.  Num.  Meth.  Fluid.<.  accepted. 

C.  PAPERS  SUBMITTED 

24.  D.  S.  Malkus.  J.  A.  Nohel.  and  B.  J.  Plohr.  Analysis  of  new  phenomena  in  shear  flow 
of  non-Newtonian  fluids,  preprint  (19S9).  Siam.  J.  Appl.  Math.,  submitted. 

25.  D.  S.  Malkus.  J.  A.  Nohel.  and  B.  J.  Plohr.  Quadratic  dynamical  systems  describing 
phenomena  in  shear  flow  of  non-Newtonian  fluids,  preprint  (19S9),  Proceedings  IMA 
Workshop  on  Xonlinear  Evolution  Equations  That  Change  Type  (19S0),  Springer 
\'erlag  Lecture  Notes,  submitted. 

2G.  J.  A.  Nohel.  R.  Pego.  and  A.  E.  Tzavaras.  Stability  of  discontinuous  steady  states 
in  shearing  motions  of  a  non-Newtonian  fluid,  preprint  (19S9),  Proc.  Royal  Soc. 
Edingburgh.  Series  A.,  submitted. 

27,  J.  A.  Nohel.  R.  Pego.  and  A.  E.  Tzavaras.  Nonlinear  stabilit}'  in  non-Newtonian  flows, 
preprint  (19S9).  Proceedings  IMA  WorLshop  on  Xonlineai  Haves  (lOSOj,  Springer 
\’erlag  Lecture  Notes,  submitted. 

D.  PAPERS  IN  PREPARATION 

2S.  D.  Brandon.  Global  e.xistence  and  asymptotic  stability  for  a  nonlinear  integrodiffer- 
ential  equation  modeling  heat  flow  (in  prepa.ation). 

29.  D.  Brandon  and  W.  J.  Hrusa.  Global  existence  of  smooth  shearing  motions  of  vicoelas- 
tic  fluids  fin  preparation). 

30.  M.  W.  Johnson  and  D.  S.  Malkus.  Shear-flow  instabilities  and  numerical  algorithms 
for  non-Newtonian  flows  (in  preparation). 

31.  D.  S.  Malkus  and  Y.-C.  Tsai.  Stability  analysis  of  implicit-explicit  time  integration 
for  viscoelastic  flow  (in  preparation). 

32.  D.  S.  Malkus  and  Y.-C.  Tsai.  New  algorithms  for  non-Newtonian  Flows  (in  prepara¬ 
tion  '). 

33.  M.  Slemrod  and  A.E.  Tzavaras.  A  viscosity  approach  for  the  .solution  of  the  Riemann 
problem  in  Eulerian  gas  dynamics  (in  preparation). 

34.  A.E.  Tzavaras.  Interplay  of  thermal  softening  and  strain  rate  sensitivity  on  the  re¬ 
sponse  of  shearing  motions  (in  preparation). 
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E,  THESES 

3o.  K.  Barki,  Ph.  D.  Thesis,  Dept  Engineering  Mechanics.  University  of  Wisconsin,  Madi¬ 
son  (in  preparation),  ab.stract  follows: 

THESIS  ABSTRACT 

The  thesis  summarizes  the  recent  development  of  viscous  material  models  and  proposes 
a  vis.-oelastic-viscoplastic  mateiial  model  which  includes  a  widely  used  viscoplastic  and 
viscoelastic  models  a^  special  cases.  The  model  also  facilitates  the  inclusion  of  a  variety  of 
internal  variables  proposed  recently  by  other  researchers.  A  modified  Euler  algorithm  is 
used  rc'  integrate  the  constitutive  equation.  An  implicit-explicit  predictor  corrector  scheme 
is  usi‘d  to  solve  the  transient  equations.  Derivation  of  stability  conditions  for  viscoelastic 
part  is  illustrated  using  the  central  difference  method.  The  development  of  the  general 
purpose  iionl'mear  program  is  discusse'-’.  The  use  of  supercomputing  also  is  discussed  in 
the  thesis,  although  the  thesis  concerns  itself  principally  with  the  material  behavior  of  the 
model.  Xumerical  e.xamples  for  one  dimensional  and  plane  strain  models  are  presented 
with  a  variety  of  loadings. 


5.  Additional  References 

.1.  Oldroyd.  Xon-Xewtonian  effects  in  steady  motion  of  some  idealized  elastico-viscous 
licpiids.  Proc.  Roy.  Soc.  London.  A  245  (195S)  278-297. 

G.  \’inogradov.  A.  Malkin,  Yu.  Yanovskii.  E.  Borisenkova.  B.  Yarlykov,  and  G.  Berezh- 
naya.  \'iscoelasfic  properties  and  How  of  narrow  distribution  poljujutadienes  and  polyiso- 
prenes  J.  Polymer  Sci..  Part  A-2.  10  (1972). 

M.  Johnson  and  D.  Segalman.  A  Model  for  viscoelastic  fluid  behavior  which  allows  non- 
afflne  deformation.  J.  Non-Xe'^'tonian  Fluid  Mech..  2  (1977)  255-270. 

J.  Pearson.  Mechanics  of  Polymer  Processing,  Elsevier  Applied  Science. London  (19S5). 

M.  Denn.  Issues  in  viscoelastic  fluid  dynamics.  Annual  Reviews  of  Fluid  Mechanics.  1989, 
to  appear. 


6.  Professional  Personnel 

D.  S.  Malkus.  Profe.ssor.  Engineering  Mcchaiaics  ^  Center  for  the  Math.  Sciences 

J.  A.  Xohel,  Professor.  Mathematics  Sv:  Center  for  the  Math.  Sciences 

R.  C.  Rogers.  Van  Vleck  Assistant  Professor.  Mathematics  <kj  Center  for  the  Math.  Sciences 

E.  Tzavaras.  Van  Vleck  Assistant  Professor,  Mathematics  Center  for  the  Math. 
Sciences 

Debor.ah  Brandon,  Postdoctoral  Researcher,  summer  1988.  currently  at  Institute  for  Math¬ 
ematics  and  Its  Applications,  University  of  Minnesota 

K.  barki,  Ph.  D.  candidate.  Engineering  Mechanics 
Y.-C.  Tsai.  Ph.  D.  candidate.  Engineering  Mechanics 
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7.  Interartions 


Tlu"  tollowiag  presentations  involving  significant  research  interactions  were  given  during 
the  repijrting  period: 

D.  S.  Malkus,  Type  change  in  non-Xewtonian  fluids.  Fluids  Research  Oriented  Group 
St'ininar.  Michigan  Technological  Institute,  Houghton,  jMI.  May,  19S7. 

D.  S.  Malkus,  Xunierical  aspects  of  the  High  Weissenberg  X'uinber  Problem,  oth  Interna¬ 
tional  worksln)p  on  Xumerical  Methods  in  X’on-X’ewtonian  Flows,  Lake  .4rrowhead.  C-4. 
June,  19S7, 

D.  S.  Malkus,  Xew  transient  algorithms  for  aon-XMwtonian  flows.  Mathematics  Department 
Seminar,  Illinois  Institute  of  Technology,  Xov..  19S7. 

R.  C.  Rogers,  “Xonlocal  materials",  MIP.4C  Seminar,  UW-Madison.  12/2/S7. 

-4,  E,  Tzavaras,  Invited  lecture  based  on  joint  work  with  Xohel  and  Rogers.  Special  Session 
on  oherra  Iiuegral  Equations",  combined  Midwest  and  Southeast  Differential  Equations 
Conference.  Vanderbilt  University,  Xashville,  TX,  Oct.  19S7. 

B,  ,1.  Plohr,  based  on  joint  work  with  D.  S.  Malkus  and  J.  A.  Nobel  ’‘Instabilities  in  shear 
flows  of  viscoelastic  fluids  with  fading  memory”,  Seminar  on  PDE's  and  Con^'inuum  Models 
of  Phase  Transitions.  Nice.  France.  1/1S-22/S8. 

D,  S.  -Malkus.  "X'ew  transient  algorithms  for  non-Xewtonian  flows".  Ill,  Inst,  of  Tech..  Math 
Department  semiiiar.  1/26/SS. 

D.  S.  Malkus.  .4  stress  calculator  for  BK2Z  and  similar  single  integral  constitutive  equa¬ 
tions.  .4nnual  Meeting  of  the  Society  of  Rheology,  .4tlanta.  Feb.,  19SS. 

J.  -4.  Xohel.  "Mathematical  problems  in  viscoelasticity".  Rheology  Research  Center,  UW- 
Madi,son,  March  4,  19SS. 

.4.  E.  Tzavaras.  Invited  lecture  b-ised  on  joint  work  with  Xohel  and  Rogers.  Special  Session 
on  "Qualitative  Theory  of  Nonlinear  PDE",  S43rd.  Meeting  of  the  AMS,  University  of 
Maryland.  College  Park.  MD,  April  19SS. 

.4.  E.  Tzavaras,  ‘‘Formation  of  shear  bands”,  LCDS  SemincU*  in  .4pplied  Mathematics, 
Brown  University,  April  19SS. 

D.  S.  Malkus.  based  on  joint  work  with  J.  A.  Xohel  and  B.  J.  Plohr  ,  ‘‘Shearing  flows  of  non- 
Xewtonian  fluids".  Si.xth  .4rmy  Conference  on  .4ppUed  Math,  and  Computing,  Boulder. 
CO.  May  31.  19SS. 

,4.  E.  Tzavaras.  Invited  participant  and  speaker,  meeting  on  ‘‘Hyperbolic  Systems  of  Con¬ 
servation  Laws",  Mathcmatisches  Forschungsinstitut.  Oberwolfach.  West  Germany,  July 
19SS. 

J  .4.  Xohel.  Seminar.  LeCchetz  Center  for  Dynamical  Systems,  Brown  University,  July  14. 
19S3. 

J.  .4.  Xohel.  Invited  lecture,  AMS  Summer  Research  Conference  on  ‘‘Current  Progress  in 
Hyperbolic  Systems:  Riemann  Problems  and  Computations",  Bowdoin  College,  Brunswick, 
Maine.  July  lG-22.  19SS. 
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J.  A.  Xohel,  Invited  lecture.  International  Conference  on  "Problems  Involving  Change  of 
Type".  University  of  Stuttgart,  October  11-14.  19SS. 

A.  E.  Tzavaras,  Seminar,  University  of  Michigan,  November,  19SS. 

J.  A.  Xohel.  .4pplied  Math.  Colloquium.  Princeton  University,  December  16,  19SS. 

,1.  .A.  Xohel,  Invited  participant  in  three  Workshops  on  nonlinear  problems  at  Institute  for 
Mathematics  and  Its  Applications.  University  of  Minnesota.  Winter-Spring,  1989. 

.4.  E.  Tzavaras.  Invited  participant,  program  on  Nonlinear  Waves  at  Institute  for  Mathe¬ 
matics  and  Its  .Applications,  University  of  Minnesota,  Winter-Spring,  1989. 

.4.  E.  Tzavaras.  Two  lectures  at  Division  for  .Applied  Mathematics,  Brown  University, 
.lanuary  1989. 

D.  S.  Malkus.  "Implicit -explicit  time  integration  methods  for  shear  flows  wdth  spurt  in¬ 
stabilities."  60th  .Annual  Meeting  of  the  Society  of  Rheology,  Gainsville,  Fla.,  Feb.  13, 
19S9. 

During  the  reporting  period  professional  interactions  (in  some  cases  joint  work)  contin¬ 
ued  with  B.  J.  Plolrr  (UW-Madison),  W.  J.  Hrusa  (CMU),  L.  Tartar  (CMU),  M.  Slemrod 
(UW-Madison).  C.  M.  Dafermos  (Brown),  R.  Pego  (U.  of  Michigan),  M.  Renardy  (VPI), 
D.  D.  Joseph  (U.  of  Minnesota),  E.  T.  Olsen  k  B.  Bernstein  (IIT),  T.  J.  R.  Hughes  (Stan¬ 
ford).  James  Walters  (BRL).  .Arthur  R.  Johnson  (Leader.  Structural  Mechanics  Group. 
Watertown),  and  members  of  the  Madison  Rheology  Research  Center. 
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